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MOLECULAR  MODELING  OF  PMN  CERAMICS 


INTRODUCTION 

Electrostrictive  materials  are  capable  of  delivering  huge  strains  when  stressed 
by  electric  fields.  One  of  the  leading  candidate  electrostrictive  materials  is  lead 
magnesium  niobate  (PMN),  which  is  currently  attracting  much  interest  in  the 
transducer  community.  The  strains  that  PMN  can  achieve  on  application  of  high 
field  are  greater  than  the  strains  in  normal  piezoelectric  ceramics  such  as  lead 
zirconate  titanate  (0.1%).  It  is  also  noteworthy  that  PMN  does  not  have  to  be  poled 
as  does  conventional  piezoelectric  ceramic.  It  belongs  to  a  class  of  relaxor 
ferroelectrics  which  displays  a  diffuse  phase  transformation.  However,  in  contrast  to 
normal  ceramics  that  are  capable  of  undergoing  a  phase  transformation  into  a  long- 
range  ferroelectric  state,  a  relaxor  such  as  PMN  displays  micropolar  regions  about  50 
A  in  diameter.  To  rationalize  the  pronounced  electrostrictive  effect  in  PMN 
materials,  Cross^  has  examined  the  role  of  compositional  heterogeneity  of  fluctuations 
in  A(BiB2)03  relaxor  perovskites  where  A,  B^  and  Bj  are  cations  (in  the  case  of  PMN, 

A  is  lead,  B^  is  magnesium,  and  Bj  is  niobium).  In  the  Cross  model,  compositional 
variations  in  B-site  ordering  lead  to  differences  in  local  and  global  symmetry.  Since 
local  symmetry  is  lower  than  the  global  symmetry,  polar  orientation  states  of  each 
microregion  are  not  identical,  and  the  energies  of  each  polar  state  are  no  longer 
degenerate.  This  leads  to  the  creation  of  dynamically  disordered  polarization  states.  The 
ramifications  of  dynamically  disordered  polarization  are  far  reaching,  for  this  model 
implies  that  coupled  dipoles  will  oscillate  and  respond  to  an  applied  electric  field 
leading  to  large  and  diffuse  permittivity  and  large  electrostrictive  strains. 

Ordering  of  the  B-site  atoms  is  in  part  governed  by  the  ionic  radii  and  charge 
of  the  B  cations;  thus,  the  tendency  to  order  is  due  to  elastic  forces  and  electrostatic 
interactions.  Using  transmission  electron  microscopy  (TEM),  Viehland^"^  has 
identified  ordered  regions  in  Pb(Mgi/3Nb2/3)03  compositions  where  the  B-site  ion 
ratio  is  1:2.  During  annealing,  these  regions  do  not  coarsen.  The  inability  of  the 
micropolar  regions  to  grow  in  size  during  annealing  is  due  to  microstructural 
chemical  effects  that  prevent  the  development  of  the  normal  long-range  order. 
However,  in  a  material  such  as  Pb(Sci/2Tai/2)^3'  which  has  an  ordered  B-site  ratio  of 
1:1,  the  ordered  domains  coarsen  during  annealing.'*  This  observation  led  Viehland 
and  others  to  suggest  that  the  clustered  regions  in  PMN  are  characterized  by  a 
nonstoichiometric  1:1  ordering  of  Mg  and  Nb.^'*  With  a  1:1  ratio  of  B-site  ions,  the 
ordered  region  is  electrically  nonneutral.  The  development  of  the  ordered  regions  is 
then  dictated  by  opposing  stabilizing  long-range  and  destabilizing  short-range 
electrostatic  interactions.  Viehland  has  suggested  that  the  size  of  the  ordered 
domains  is  a  result  of  charge  variations  due  to  mixed  B-site  occupancy.  According  to 
this  view,  ordered  domain  size  and  subsequently  the  electromechanical  properties  can 
be  manipulated  by  creating  defect  structures  such  as  can  be  achieved  by  the 
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introduction  of  oxygen-deficient  sites.  Substitutions  for  A-site  atoms  can  bring  about 
changes  in  annealing  and  ordering  behavior.  For  example,  substitution  of  La^  for 
Pb'^  forces  ordering  of  Nb^""  and  However,  substitution  of  Ba^""  or  for  Pb^^ 

significantly  reduces  the  relaxor  response.*  It  is  clear  that  the  A-site  ion  must 
influence  the  size  of  the  ordered  domains  and  therefore  the  electromechanical 
response.  Because  La*"^  is  aliovalent,  its  substitution  into  the  A-site  may  be 
accompanied  by  the  formation  of  oxygen  vacancies  to  compensate  for  the  positive 
charge. 

Another  complication  in  the  PMN  puzzle  has  to  do  with  the  possibility  of 
structural  nonuniformities  in  the  ordered  1:1  microregion.  Experiments  employing 
atomic  resolution  transmission  electron  microscopy  have  demonstrated  the  existence 
of  strong  nonuniform  polarized  transverse  strain  waves  in  the  ordered  regions.’  One 
hypothesis  for  this  observation  is  that  the  nonuniform  distortions  are  due  to  the 
presence  of  an  "embryonic"  Pb2Nb207  pyrochlore  phase  with  the  perovskite  phase. 

The  presence  of  this  inhomogeneity  may  be  responsible  for  the  development  of  local 
polarization  and  the  absence  of  a  long-range  ferroelectric  phase. 

Quite  germane  to  PMN  is  an  alternative  atomistic  model.  Thomas  proposed 
that  indirect  <110>  coupling  dominates  in  PMN  and  influences  B-site  ordering.*  In 
the  Thomas  model,  the  temperature  at  which  an  individual  octahedron  becomes 
decoupled  from  its  neighboring  octahedron  gives  rise  to  a  Curie  point  distribution 
function  which  can  explain  the  diffuseness  of  the  ferroelectric-paraelectric  phase 
transitions  in  PMN  ceramics.  A  ferroelectrically  active  microregion  consisting  of 
active  NbO^  octahedra  is  surrounded  by  inactive  MgO^  octahedra.  Between  the  NbO^ 
cluster  and  MgO^  surroundings  there  exists  a  dynamic  wall  rich  in  MgOg.  The  NbOg 
octahedra  within  the  wall  are  coordinated  (coupled)  with  fewer  NbO^  so  that  they 
undergo  a  ferroelectric-paraelectric  transition  and  decoupling  at  a  lower  temperature 
than  the  coupled  NbO^  within  the  cluster.  As  decoupling  within  the  wall  proceeds, 
more  NbOg  octahedra  become  uncoupled  and  the  "wall"  propagates  towards  the 
center  of  the  microregion.  All  of  this  can  give  rise  to  the  diffuse  phase  transition 
characteristic  of  PMN  ceramics. 

It  is  therefore  abundantly  clear  that  to  manipulate  the  el^tromechanical 
response  of  PMN,  attention  must  be  given  to  compositional  control  of  PMN. 
Compositional  modification  can  be  achieved  by  altering  the  A-  or  B-site  atom  or  by 
oxygen-deficient  ordering.  These  modifications  can  lead  to  optimization  of  the 
electromechanical  properties  of  PMN.  Current  research  has  been  directed  to  material 
development  of  PMN  compositions  with  special  emphasis  on  improved  performance 
such  as  increased  induced  strain,  decreased  strain  hysteresis,  increased  coupling,  and 
reduced  temperature  dependence  of  strain  and  polarization  behavior. 

The  objective  of  the  present  research  was  to  examine  the  effects  of  A-site 
substitutionals  on  the  bonding  and  stability  of  the  ordered  regions,  from  the 
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perspective  of  semiempirical  quantum  mechanics.  We  propose  that  ordering 
behavior  in  PMN  might  be  related  to  interactions  between  the  orbitals  of  the  A-site 
atom  with  the  orbitals  involved  in  bonding  of  the  B-site  atoms  and  oxygen.  Because 
there  has  been  no  analysis  of  orbital  interactions  in  PMN  ceramics,  it  may  be  fruitful 
to  examine  how  orbital  overlap  or  bonding  between  Nb  and  O  and  between  Mg  and 
O  is  related  to  the  nature  of  the  A-site  atoms.  There  are  several  ways  to  depict  these 
orbital  interactions.  In  a  NbO^'  perovskite,  niobium  is  coordinated  in  the  <100> 
direction  in  an  octahedral  environment  by  six  oxygen  atoms.  Niobium-oxygen 
bonding  can  involve  the  interaction  of  2s  and  2p  orbitals  in  oxygen  and  the  eg 
and  4dx2.y2)  orbitals  in  niobium  (figure  1).  The  remaining  niobium  tjg  (4d^,  and 
4dy^ )  orbitals  might  be  expected  to  interact  with  the  2p  oxygen  orbitals  to  generate  k- 
bonding  (figure  1).  We  are  hypothesi2dng  that  the  nature  of  the  A-site  atom 
irdluences  the  local  bonding  in  the  PMN  crystal  lattice,  and  this  is  brought  about  by 
changing  the  electron  count  and  position  of  the  Fermi  level.  If  indeed  the  stability  of 
the  PMN  microregions  is  related  to  bonding,  can  a  predictive  model  be  formulated  to 
show  how  substitution  of  other  atoms  on  the  A-site  is  related  to  the  size  and 
character  of  the  microregions  necessary  for  electrostriction? 

This  proposal  represents  an  effort  to  explain  the  electrostrictive  phenomenon 
in  PMN  ceramics  within  the  framework  of  molecular  orbital  theory,  which  has  been 
used  by  chemists  to  successfully  predict  and  study  structure-property  relationships  of 
numerous  solid-state  materials.  In  effect,  the  research  suggested  here  will  bring 
together  the  language  of  chemistry  and  physics  to  study  relaxor  properties.^^’^^'^  An 
atomistic  approach  to  the  composition  variation  in  PMN  and  the  implications  of 
composition  in  electrostriction  may  ultimately  result  in  a  predictive  model.  The  work 
described  here  is  based  on  a  previous  study  of  the  effects  of  orbital  interactions  in 
ordered-defect  rock-salt  structures,  where  it  was  demonstrated  that  metal-metal  and 
metal-oxide  bonding  helps  to  stabilize  these  structures.^*  The  approach  in  the  present 
study  was  to  determine  the  role  of  band  structures,  density  of  states,  and  orbital 
overlap  populations  on  bond  stability  in  model  perovskite  structures  representative  of 
PMN. 


THEORY 


An  understanding  of  orbital  interactions  in  solid  materials  can  lead  to  an 
improved  understanding  of  material  properties  such  as  bond  stability.  Bonding  in 
crystal  structures  can  be  formulated  in  the  language  of  band  orbitals.  The  theory  is 
described  in  the  next  section. 
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BAND  ORBITALS 


Orbital  interactions  in  several  crystal  models  of  PMN  were  studied  using 
extended  Hiickel  (EH)  molecular  orbital  theory.  When  applied  to  solid-state 
materials,  the  EH  approach  gives  a  band  structure  calculation,  which  can  be  used  to 
analy2!e  bonding.  When  applied  to  a  one-dimensional  lattice  containing  one  atom  per 
unit  cell,  the  atomic  orbital  (cp)  can  be  represented  by  a  Bloch  orbital  shown  as 
follows: 

=N exp(j^)X;i,  (1) 

n 


where  is  an  atomic  orbital  on  the  atom  in  the  n-th  unit  cell,  N  is  the  number  of 
unit  cells  in  the  solid,  k  is  a  wavevector  along  the  one-dimensional  direction,  and  a  is 
the  lattice  dimension  of  the  one-dimensional  cell.  Values  of  the  wavevector  are 
defined  by 


(2) 


Crystal  orbitals  are  linear  combinations  of  Bloch  orbitals: 

(1=1 

The  eigenvalues  of  the  Hamiltonian,  H,  of  equation  (3)  represent  the  energies,  Ej  (k), 
of  the  crystal  orbitals: 


iP¥,ik)=E,m^k). 


Minimizing  the  energy  with  respect  to  the  coefficients  gives  a  secular  determinant; 


where 
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and 


(7) 


Solution  of  equation  (5)  leads  to 


E(k)=a+2^coska, 

where  a  is  the  Coulomb  integral  and  P  the  nearest  neighbor  resonance  integral.  For 
a  series  of  overlapping  orbitals,  a  plot  of  E(k)  vs.  k  can  be  obtained.  The  width  of 
the  band,  or  the  magnitude  of  interaction  between  the  orbitals  on  neighboring  atoms, 
is  given  by  4p.  In  other  words,  a  chain  of  strongly  interacting  orbitals  results  in  a 
large  dispersion  or  large  4p. 

The  plot  of  E(k)  vs.  k,  which  is  called  a  band  structure  diagram,  allows  one  to 
analyze  the  electronic  properties  of  solid-state  materials.  The  range  of  k  is  called  the 
first  Brillouin  zone,  which  can  be  constructed  using  the  reciprocal  crystal  lattice. 
From  the  band  structure,  we  can  calculate  the  density  of  states  (DOS),  which 
represents  the  number  of  levels  at  a  given  energy.  To  calculate  the  bonding 
characteristics  between  two  atoms  in  a  solid-state  material,  the  overlap  population 
represents  the  overlap  integrals  summed  over  all  orbitals  of  the  two  atoms.  The 
overlap  population  weighted  density  of  states  is  called  the  crystal  orbital  overlap 
population  (COOP).  The  COOP  describes  the  way  the  orbitals  interact  and  allows 
one  to  determine  the  strength  of  bonding. 


MODEL  CRYSTALS 

The  crystal  structure  of  PMN  at  room  temperature  is  a  cubic  perovskite 
belonging  to  space  group  Ptn3m.  The  lattice  dimension  is  about  4.04  To  study 

the  effects  of  A-site  substitution  on  Nb-O  bonding  along  the  <100>  crystallographic 
directions,  the  models  chosen  for  this  study  were  PbNbOj^"^  and  LaN^a^"^  .  Figure  2 
shows  the  unit  cell  crystal  structure  of  PbNbOj^*.  The  structure  of  LaNbOa^  is 
similar.  These  structures  have  Nb  positioned  at  the  corners  (0,  0,  0)  of  a  primitive 
cubic  unit  cell.  The  Nb  is  coordinated  along  the  cell  edges  with  oxygen  atoms  at  the 
edge  centers  at  (1/2,  0,  0;  0,  1/2,  0;  0,  0,  1/2).  The  A-site  cation,  Pb  or  La^,  is 
placed  in  the  body-centered  position  at  (1/2,  1/2,  1/2).  The  size  of  the  PbNl^a^* 
and  LaNbOa'^  unit  cell  was  set  at  3.7  A,  which  implies  a  bond  length  of  1.85  A  for 
Pb-O  and  La-O.  This  size  was  chosen  to  correct  for  the  diffuseness  of  the  atomic 
orbitals  and  to  maximize  orbital  overlap. 
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The  crystals  were  modeled  as  three-dimensional  lattice  systems  with  the  k 
vector  changing  within  a  three-dimensional  Brillouin  zone.  The  points  in  the  three- 
dimensional  zone,  which  correspond  to  certain  points  of  the  cubic  lattice  as  shown  in 
figure  3,  are  T,  M,  X,  and  The  directions  in  reciprocal  space  in  this  cubic  lattice 
are  therefore  T  to  M,  M  to  X,  X  to  R,  and  R  to  T.  k-space  meshes  consisted  of  40  k 
points.  The  Huckel  parameters  used  in  the  molecular  orbital  calculations  are  listed  in 
table  1. 


RESULTS  AND  DISCUSSION 


A  preliminary  study  was  first  undertaken  to  obtain  a  qualitative  orbital 
description  of  bonding  in  the  model  crystal  structures.  The  electronic  configurations 
of  Nb,  Pb,  and  La  are  listed  in  table  2.  Niobium  uses  its  4d  and  5s  electrons  to  form 
a-bonds  and  7t-bonds  with  the  2s  and  2p  orbitals  of  surrounding  oxygen  atoms.  To 
examine  the  nature  of  this  bonding  and  how  it  changes  along  the  Brillouin  zone  and 
affects  the  band  structure,  it  is  best  to  consider  the  interactions  to  form  n-bonds  in  a 
simple  one-dimensional  repeating  chain.  Figure  4  shows  the  two  band  structures  of  a 
one-dimensional  chain  of  alternating  niobium  and  oxygen  and  the  interactions 
between  the  oxygen  2pj,  and  niobium  4dx2  Bloch  orbitals.  At  k  =  0,  the  2p^  and  4dxz 
Bloch  orbitals  do  not  mix  because  of  differences  in  symmetry.  The  2p,  Bloch  orbital 
is  symmetric  with  respect  to  rotation  about  the  z-axis,  whereas  the  4dxj  orbital  is 
asymmetric.  This  implies  that  these  orbitals  should  not  interact  at  k  =  0.  However, 
at  k  =  7r/2,  both  orbitals  are  allowed  to  mix  because  of  symmetry.  As  a  result  of 
orbital  mixing,  two  new  states  are  created,  one  characterized  by  strong  orbital 
overlap  and  therefore  at  lower  energy  and  the  other  having  antibonding  interactions. 
The  direct  consequence  on  changing  the  wavevector  from  T  to  X  leads  to  a  gradual 
lowering  of  the  energy  of  the  2p2  Bloch  orbital  and  a  gradual  increase  in  4dx2. 

Similar  arguments  can  be  invoked  for  the  more  complex  three-dimensional  structures. 
Although  the  resulting  band  structures  can  be  quite  complicated,  orbital  interaction 
can  be  depicted  as  in  the  case  of  the  simple  one-dimensional  repeating  chain  to  gain 
an  intuitive  perception  of  orbital  interactions.  The  number  of  bands  in  the  structure 
is  equivalent  to  the  number  of  interacting  orbitals.  Some  bands  will  be  steep  because 
of  strong  orbital  overlap  and  others  will  be  flat  The  bands  are  influenced  by  both 
the  symmetry  and  topology  of  the  interacting  orbitals.  In  general,  more  diffuse 
orbitals  such  as  s  and  p  give  rise  to  strong  dispersion;  contracted  orbitals  such  as  d 
result  in  flat  and  dispersionless  bands. 

It  is  more  instructive  to  analyze  DOS  and  COOP  plots  rather  than  the  complex, 
sphaghetti-like  band  structures.  The  DOS  plots  of  PbNbOj^*,  and  LaNbOj^  are 
illustrated  in  figures  5  and  6,  respectively.  The  separate  orbital  contributions  of  the 
oxygen,  A-site  lead  or  lanthanum,  and  B-site  niobium  to  each  DOS  are  shown  as 
dashed  lines  in  figures  7  to  21.  Consider  first  the  DOS  plots  of  PbNb03  .  The 
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interactions  between  the  5p,  6s,  and  eg  and  4d^2_y2  )  orbitals  of  niobium  with  the 

2p  orbitals  of  oxygen  contribute  to  a-  and  Tt-bonding  (figure  1).  This  is  shown  by 
the  contributions  to  the  overall  DOS  of  the  O  2p  bands  in  figure  7  and  of  the  Nb  5s, 
5p,  and  eg  bands  in  figures  8  to  11.  The  formation  of  7t-bonding  can  occur  by 
overlap  between  the  5p  and  t2g  niobium  orbitals  and  2p  oxygen  orbitals.  This  is 
illustrated  in  figure  1,  which  shows  the  overlap  between  a  d^^  niobium  and  a  2p 
oxygen  orbital.  Similar  considerations  apply  to  the  bands  in  LaNbOj^*  (figures  12  to 
16). 


Now,  let  us  consider  the  contributions  of  the  A-site  atomic  orbitals  to  the  DOS. 
For  PbNb03^*,  the  contributions  of  the  6s  and  6p  lead  orbitals  to  the  DOS  are  shown 
in  figures  17  and  18.  Here  the  6s  and  6p  bands,  which  are  well  separated  and  at 
higher  energies,  make  a  small  contribution  to  the  density  of  the  2p  oxygen  and  5s 
and  5p  niobium.  The  contributions  of  the  lanthanum  6s,  6p,  and  tjg  orbitals  are 
exemplified  in  figures  17  to  21.  Again,  we  note  small  contributions  to  the  lower 
energy  bands.  We  conclude  that  A-site  atoms  do  not  bond  with  neighboring  atoms 
in  the  crystal  lattice. 

So  far,  our  analysis  has  revealed  nothing  unusual  from  the  DOS  plots.  The 
orbital  interactions  to  form  a-  and  ;i-bonds  support  our  expectations.  The  nature  of 
the  A-site  atom  does  not  seem  to  influence  these  interactions.  For  a  complete  picture, 
however,  we  need  a  way  to  describe  the  bands  as  either  bonding  or  antibonding.  We 
can  do  this  by  calculating  the  overlap  integrals  between  the  niobium  and  oxygen 
summed  over  all  orbitals  in  these  atoms  to  give  a  COOP.  Positive  overlaps  signify 
bonding;  negative  overlaps  imply  antibonding.  As  shown  in  figures  22  and  23,  the 
COOP'S  of  PbNbOj^*  and  LaNbOj’*  are  not  identical.  To  appreciate  this  difference, 
we  have  specified  the  position  of  the  Fermi  level  for  both  models.  The  Fermi  level 
lies  midway  between  the  highest  bonding  and  lowest  antibonding  orbitals.  Its 
position  is,  therefore,  determined  by  the  number  of  electrons  occupying  the  orbitals. 
The  trivalent  state  of  lanthanum  has  one  electron  less  than  the  divalent  state  of  lead. 
This  suggests  that  the  Fermi  level  should  be  lower  when  La^*  is  substituted  in  the  A- 
site.  Indeed,  calculations  show  that  in  LaNbO,**,  the  Fermi  energy  lies  at  about  -14.5 
eV,  whereas  in  PbNbOa^*  the  Fermi  level  is  -10.1  eV.  It  is  particularly  significant  that 
near  the  Fermi  level  in  PbNbOj^'’  there  is  an  antibonding  Nb-O  band  partially 
occupied  by  electrons.  This  suggests  that  substitution  of  Pb^*  in  the  A-site  may 
destabilize  the  bonding  between  niobium  and  oxygen  giving  rise  to  weaker  bonds  in 
the  crystal  structure.  The  magnitude  of  this  destabilization  is  such  that  long-range 
ordering  and  crystallization  may  be  more  difficult  when  Pb^*  exclusively  occupies  the 
A-site.  Substitution  of  La^"^  in  the  A-site,  however,  is  predicted  to  favor  bonding 
along  the  <100>  crystallographic  directions  and,  therefore,  ordering  and 
crystallization. 

The  calculations  in  this  investigation  must  be  regarded  as  tentative.  More 
work  is  needed  to  test  this  theory  on  various  models  and  establish  trends.  To 
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examine  the  consequences  of  bonding  and  orbital  interactions  in  PMN,  it  is  proposed 
to  extend  this  preliminary  work  to  superlattice  structures  of  PMN,  to  examine  the 
orbital  and  electronic  effects  of  substitutions  of  Ba-^  Sr^^  and  other  cations  in  the  A- 
site,  and  to  correlate  these  effects  with  the  position  of  the  Fermi  level. 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  original  objective  of  this  work  was  to  demonstrate  that  molecular 
modeling  is  a  practical  tool  for  the  study  of  ceramics  used  in  Navy  underwater 
transducers.  Given  the  formidable  complexity  of  the  ceramic  PMN,  there  is  strong 
reason  to  develop  structure/ property  relationships,  which  may  lead  to  a  judicious 
selection  of  PMN  compositions  and  processing  conditions.  We  chose  to  study  the 
effect  of  an  A-site  substitutional,  lanthanum,  on  the  annealing  properties  of  PMN,  a 
candidate  electrostrictive  ceramic  for  high  power  sonar  projection.  To  approach  this 
problem,  we  have  chosen  a  semiquantum  mechanical  approach  called  extended 
Hiickel  theory  where  band  structures,  density  of  states,  and  crystal  overlap 
populations  can  be  calculated.  These  calculations  can  be  used  to  explain  certain 
properties  of  the  ceramic.  In  the  present  work,  we  propose  a  tentative  explanation 
for  the  lack  of  ordering  in  PMN  ceramic  when  lead  exclusively  occupies  the  A-site  in 
the  crystal  lattice.  In  contrast,  when  lanthanum  is  substituted  for  lead,  there  is  an 
enhancement  in  bond  stability  and  a  tendency  for  long-range  ordering  and  annealing. 
This  markedly  changes  the  electrostrictive  properties  of  the  ceramic. 

Of  course,  this  is  only  part  of  the  story.  We  have  not  considered  the  potential 
role  of  oxygen  defects  in  La-substituted  PMN,  the  effects  of  other  substitutionals,  and 
the  nature  of  the  electrostatic  effects  arising  because  of  charge  imbalances  between 
the  ordered  and  nonordered  regions.  These  are  all  issues  which  will  have  to  undergo 
scrutiny.  In  the  meantime,  we  are  extending  this  work  to  develop  a  more  complete 
picture  of  the  role  of  orbital  interactions  in  ceramics  for  transducers. 
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Figure  1.  Orbital  Interactions  in  PMN  Ceramics 


Figure  2.  Unit  Cell  Crystal  Structure  of  PbNbOj* 
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Figure  4.  Band  Energies  of  One-DimensioUal  Chain  of  Niobium  and  Oxygen  Atoms 
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Figure  6.  Density  of  States  (DOS)  of  LaNbO; 
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Figure  7.  Contribution  of  Oxygen  2p  Bands  to  DOS  in  PbNbO- 
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5p  Bands  to  DOS  in  PbNbO, 


m 


m 


lO 


o 

CM 

I 


A9  'AO>I3N3 


18 


A3  'AoaaNH 


19 


U5 


lO 


lO 


O 

CM 

I 


A9  'AO>I3N3 


20 


Figure  13.  Contribution  of  Niobium  5s  Bands  to  DOS  in  LaNbOj 
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Figure  14.  Contribution  of  Niobium  5p  Bands  to  DOS  in  LaNbOj 
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Bands  to  DOS  in  LaNbO^ 
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Figure  16.  Contribution  of 
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6s  Bands  to  DOS  in  PbNbO. 
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Figure  19.  Contribution  of  Lanthanum  6s  Bands  to  DOS  in  LaNb03 
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Figure  21.  Contribution  of  Lanthanum  tig  Bands  to  DOS  in  LaNbO: 
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Figure  23.  Crystal  Orbital  Overlap  Population  (COOP)  in  LaNbO: 


Table  1.  Extended  Huckel  Parameters® 


atom 

orbital 

Hu 

eV 

B 

C2’’ 

B 

C2 

Nb 

5s 

-10.1 

1.89 

- 

1.0000 

- 

5p 

-6.86 

1.85 

- 

1.0000 

- 

4d 

-12.1 

4.08 

1.64 

0.6401 

0.55160 

Pb 

6s 

-15.7 

2.35 

- 

1.0000 

-  - 

6p 

1 

00 

o 

2.06 

- 

1.0000 

- 

La 

6s 

-7.67 

2.14 

- 

1.0000 

- 

6p 

-5.01 

2.08 

- 

1.0000 

5d 

-8.21 

3.78 

1.381 

0.77651 

0.45861 

O 

2s 

-32.3 

2.275 

- 

1.0000 

- 

2p 

-14.8 

1275 

- 

1.0000 

- 

Whangbo,  M.  Evain,  T.  Hughbanks,  M.  Kertesz,  S.  Wijeyesekera,  C.  Wilker,  C. 
Zheng,  and  R.  Hoffmann,  "Extended  Htickei  Molecular,  Crystal  and  Properties 
Package,"  QCPE  571  Program,  Quantum  Chemistry  Program  Exchange,  Indiana 
University. 

'Slater  exponents. 
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Table  2.  Electronic  Configurations^ 


atom 

electronic  configuration 

Nb 

4(P5s' 

Pb 

5d'°6s26p2 

La 

5d^6s^ 

O 

2s‘2p^ 

“Outer  orbitals  only. 
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